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We apply bulk holographic dark energy in general braneworld models with a Gauss-Bonnet 
term in the bulk and an induced gravity term and a perfect fluid on the brane. Without 
i making any additional assumptions we extract the Fricdmann equation on the physical brane 
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INTRODUCTION 



and we show that a p-p\ coupling arises naturally by the full 5D dynamics. The low- 
energy (late-time) evolution reveals that the effective 4D holographic dark energy behaves as 
"quintom", that is it crosses the phantom divide w = — 1 during the evolution. In particular, 

o ! 

the Gauss-Bonnet contribution decreases the present value of u>a, while it increases the 
?-h , growing rate of w\(z) with z, in comparison with the case where such a term is absent. 

PACS numbers: 95.36,+x, 98.80.-k, 04.50.-h 
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Holographic dark energy [1, 2, 3, 4, 5, 6] is an interesting and simple idea of explaining the 
observed Universe acceleration [7]. It arises when the more fundamental holographic principle 
[8, 9] is applied in the cosmological framework [fO, 11] (although there are some objections 
^ ■ on this approach [12]). Holographic dark energy reveals the dynamical nature of the vacuum 

h : 

energy by relating it to cosmological volumes. The background on which it is based, is the 
black hole thermodynamics [13, 14] and the connection between the UV cut-of of a quantum 
field theory, which is related to vacuum energy, and a suitable large distance of the theory [15]. 
This connection, which was also known from AdS/CFT correspondence, proves to be necessary 
for the applicability of quantum field theory in large distances. The reason is that while the 
entropy of a system is proportional to its volume the black hole entropy is proportional to its 
area. Therefore, the total energy of a system should not exceed the mass of a black hole of the 
same size, since in this case the system would collapse to a black hole violating the second law 
of thermodynamics. When this concept is applied to the Universe, the corresponding vacuum 
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energy is the holographic dark energy. 

Until now, almost all works on the subject have been formulated in the standard 4D frame- 
work. On the other hand, brane cosmology [16, 17] exhibits many phenomenological successes 
[18]. In a recent work [19] we presented a generalized and restored holographic dark energy in 
the braneworld context. The basic argument was that in such a framework black holes will in 
general be D-dimensional [13, 14] and therefore holographic dark energy should be considered in 
the bulk. Subsequently, it gives rise to an effective 4D dark energy with "inherited" holographic 
nature, and this one is present in the (also arisen from the full dynamics) Friedmann equation 
of the brane. In [19] we applied this bulk holographic dark energy in a general single-brane 
model and we reproduced the results of conventional 4D calculations [1, 2, 3, 4, 5, 6], having in 
mind that the physical interpretation is different. In [20] we applied it in a general two-brane 
model with moving branes and we showed that "quintom" behavior [21, 22, 23] arises naturally 
for a large parameter space area of a simple solution subclass, without the inclusion of special 
fields or potential terms. In particular we found that wa was larger than —1 in the past while 
its present value is w\ = —1.08, and the phantom divide wa = — 1 was crossed at z p « 0.49, a 
result in remarkable agreement with observations [24, 25]. 

In this work we examine general single-brane models, including a Gauss-Bonnet term in the 
bulk [26, 27, 28, 29, 30, 31] (see also [32, 33] for a Gauss-Bonnet term in conventional 4D 
cosmology). Such a higher-curvature combination corresponds to the leading order quantum 
correction to gravity, in an effective action approach to string theory and in particular in the 
case of the heterotic string [34] , and its coupling is related to the Regge slope parameter on string 
scale. Furthermore, the Gauss-Bonnet combination is the only curvature squared form which 
gives ghost-free self-interactions for the graviton (around flat spacetime) [35] and maintains its 
zero modes of the perturbations localized on the brane [36] . Fortunately, holographic description 
holds for braneworld Gauss-Bonnet gravity, although the subject is not trivial since there are 
some ambiguities in the case of non-flat branes away from the bulk boundary [37]. Applying bulk 
holographic dark energy in this framework, and without any additional assumption, we acquire 
the interesting situation of an interaction between the 4D dark energy and the matter density of 
the brane. In this case, cosmological evolution and in particular the dependence of the 4D dark 
energy on the brane scale factor, acquires a correction in terms of the Gauss-Bonnet coupling. 
The rest of the text is organized as follows: In section II we present the holographic dark energy 



3 



in the bulk and in section III we apply it to a general single-brane model in 4+1 dimensions 
with a Gauss-Bonnet term in the bulk. Finally, in IV we discuss the physical implications of 
our analysis and we summarize the obtained results. 



II. FORMULATION OF HOLOGRAPHIC DARK ENERGY IN A GENERAL BULK 

In this section we display the basic results of bulk holographic dark energy, formulated in 
[19]. The mass Mbh of a spherical and uncharged D-dimensional black hole is related to its 
Schwarzschild radius r s through [14, 38]: 



M BH = r?- 3 (V^M D ) D - 3 M D (2.1) 

81 K—2~) 

where the D-dimensional Planck mass Mo is related to the D-dimensional gravitational constant 
Gd and the usual 4-dimensional Planck mass M p through: 

M 2 = M£- 2 Vd_ 4 , (2.2) 

with Vd-4 the volume of the extra-dimensional space [14]. 

If Pad is the bulk vacuum energy, then application of holographic dark energy in the bulk 
gives: 

PA D Vol(S D ~ 2 ) < r D ~\^M D ) D -*M D il, (2.3) 

81 \—2-) 

where Vo\(S D ~ 2 ) is the volume of the maximal hypersphere in a .D-dimensional spacetime, given 
from: 

Vol(S D ~ 2 ) = A D r D ~\ (2.4) 

with 

D-l 
7T 2 

A D 



( D-2 \\ 



-2^7^, (2.5) 



(D-l)! 

for D — l being even or odd respectively. Therefore, by saturating inequality (2.3) introducing L 
as a suitable large distance (IR cut-off) and c 2 as a numerical factor, the corresponding vacuum 
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energy is, as usual, viewed as holographic dark energy: 

Pad = c 2 (V^M D ) D - 3 M D A^ ^ ~ 2 LT 2 . (2.6) 

As was mentioned in [19], the "suitable large distance" which is used in the definition of L in 
(2.6) could be the Hubble radius [39], proportional to the square root of the Hubble radius [4], 
the particle horizon [10], the future event horizon [1, 3, 40], or the radius of the event horizon 
measured on the sphere of the horizon [5] (see also [6] for the corresponding formulation in 
Chaplygin gas and tachyon holographic models). For a flat Universe the future event horizon 
is the most suitable ansatz and furthermore it is the only one that fits holographic statistical 
physics, namely the exclusion of those degrees of freedom of a system that will never be observed 
by the effective field theory [41]. 



III. HOLOGRAPHIC DARK ENERGY IN GENERAL 5D BRANEWORLD MODELS 
WITH A GAUSS-BONNET TERM IN THE BULK 

We are interested in applying bulk holographic dark energy in general 5D braneworld models 
with a Gauss-Bonnet term in the bulk. We consider an action of the form [26, 27]: 

S = j d^xdy^g {M\R - p A5 + MlaC GB ) + J d 4 x^ (C^ at - V + r c Mfi? 4 ) • (3.7) 

In the first integral M5 is the 5D Planck mass, /?A5 is the bulk cosmological constant which is 
identified as the bulk holographic dark energy, and R is the curvature scalar of the 5D bulk 
spacetime with metric gAB- As usual, 

C GB = R 2 - 4,R AB R AB + RabcdR ABCD (3.8) 

is the Gauss-Bonnet term with coupling constant a, and Rabcd, Rab are respectively the 
Riemann and Ricci tensors. In the second integral 7 is the determinant of the induced 4D metric 
7 Q/ 3 on the brane, V is the brane tension and 0^ at is an arbitrary brane matter content. Lastly, 
we have allowed for an induced gravity term on the brane, arising from radiative corrections, 
with r c its characteristic length scale and R4 the 4D curvature scalar [31, 42, 43]. 

In order to acquire the cosmological evolution on the brane we use the Gaussian normal 
coordinates with the following metric form [44, 45]: 

ds 2 = -m 2 (T,y)dT 2 + a 2 (T,y)dn 2 k + dy 2 . (3.9) 
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The brane is located at y = 0, we impose a ^-symmetry around it, m(r, y = 0) = 1 and 
dQi stands for the metric in a maximally symmetric 3-dimensional space with k = — 1, 0, +1 
parametrizing its spacial curvature. Although we could assume a general matter-field content 
[46], we consider a brane-Universe containing a perfect fluid with equation of state p = wp. In 
this case, and after integration of the 00 and ii components of the 5D Einstein equations around 
the brane, the low-energy (p <C V) brane cosmological evolution is governed by the following 
equation [26, 31, 42] (see also [28] for similar brane solutions): 



H 2 + ^ = (72Mf - 16a p A5 M| + 6r c VMf ) Vp + - 1 (2 + Vl+T) * , 

(3.10) 

where 

A = 2ap A5 /3Ml (3.11) 

In (3.10) a stands as usual for the brane scale factor. In order to acquire a form consistent with 
conventional 4D Friedmann equation we make the identification: 

- - (3.12) 

871-Mf or c 

and we define 

Vi(<*, PAB ) = 2 V§ Ms ^— -, (3.13) 

-2a PA5 ±M?(^-6r c ) 

where M p is the 4D Planck mass. In this case brane evolution equation (3.10) becomes: 

9 k 8tt . 8tt . . 

H + - 2 = ^ p + V^a, PA5 ) P + ^PA, (3.14) 

where the (effective in this higher-dimensional model) 4D dark energy is: 

3M„ 2 Ml t , ^\ / , ^2 

PA = PAi = — Ml-Vl + A 2 + Vl + A . (3.15 

2vr(gL^-2r c )2 96vra V M J 

In the equations above /?A5 is the 5D bulk holographic dark energy, which according to (2.6) is 
given by: 

p A5 = c 2 |-Mf L~ 2 . (3.16) 
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Relations (3.11)-(3.16) describe the low-energy (late-time) cosmological evolution on the brane. 
Similarly to [19, 20] the holographic nature of pas is the cause of the holographic nature of 
p\. Finally, the 5D Planck mass M5 is related to the standard 4D M p through M| = M 2 /L^ 
(according to (2.2)), with L5 the volume (size) of the extra dimension. 

Let us make some comments here. The above expressions in the limit a — > (where A — ► 
and Vi(a,pA5) - ► 0) tend smoothly to those analyzed in [19]. However, in the presence of the 
Gauss-Bonnet term (a 7^ 0) we observe an interesting interacting behavior. Indeed, in (3.14) 
there is a coupling between p and V\{a, pas), that is a term depending on p\5 and therefore 
on pa- We mention that the coupling between p and p\ arises naturally through the full 5D 
dynamics and the use of bulk holographic dark energy, and it is not a result of an arbitrary 
introduction by hand, which is the usual case in interacting holographic dark energy in the 
literature [39, 47] even in the case where a Gauss-Bonnet term is present [48]. 

Our final goal is to find the relation between P a and the metric scale factor a of the brane. 
However, the complex form of the above equations makes it impossible to acquire such an 
expression analytically. Therefore, in the following we describe the necessary approximations. 
Firstly, as we have already mentioned, according to (2.2) M| = M 2 /L^ with L5 the volume 
of the extra dimension. In this work we assume that L5 is arbitrary large (but not infinite), 
i.e. it is larger than any other length of the model, thus leaving brane evolution unaffected 
by the bulk size or bulk boundaries and this is the reason for the single-brane consideration. 
Therefore, in the calculations below we impose Mil Ml = L 5 » r c and 1/L 5 -> 0. The role of 
the bulk size was investigated in [20]. Secondly, we expand (3.13) and (3.15) in terms of the 
Gauss-Bonnet coupling a and we keep only the linear term. Actually this is also a consistency 
requirement since, in heterotic string theory background, the Gauss-Bonnet form is the leading 
order quantum correction to gravity, i.e we have already kept only linear terms in a [49]. These 
steps lead to: 

Vi(a, PAS ) = V!(a,L) = ^aL~ 2 + 0(a 2 ), (3.17) 

PA = 3c2 J^2 M P L ~ 2 ( 1 + a ^ L ' 2 ) + ( 3 - 18 ) 

Finally, we have to determine the cosmological length L which is present in the bulk holographic 
dark energy expression (3.16) and has been transferred to relations (3. 17), (3. 18), too. In the 
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following we will consider a flat Universe, in order to safely use the future event horizon to define 
L, without entering into the relevant discussion of the literature concerning the IR cut-off in 
non-flat cases [1, 3, 4, 5, 40]. However, the model of the present work, such as the majority 
of braneworld models of the literature, is not maximally isotropic and this feature makes the 
analytical calculation of the future event horizon an impossible task. In this anisotropic case we 
can alteratively use the 4D future event horizon (the 4D spacetime is the maximally isotropic 
subspace of the model), without losing the qualitative behavior of the observables. Fortunately, 
the calculations in the simple case without a Gauss-Bonnet term [19], showed that the use of 
the 4D future event horizon leads to identical quantitative results comparing to those obtained 
within the traditional holographic dark energy [1, 2, 3, 4, 5, 6]. 

Using the above approximations we obtain the following form for the effective 4D holographic 
dark energy: 

PA = ^Y^2 M P R h 2 (l + « ^ Rf) . (3-19) 
and substitution to Friedmann equation (3.14), for the flat-Universe case, gives: 

In these relations, the 4D future event horizon is given as usual by: 

f°° da 1 

R - = a L hp- (3 ' 21) 

Finally, we have to insert in (3.20) the known form for p(a), namely p = poa~ 3 , with po its 
present value. 

The aforementioned integral equations determine completely the brane evolution, in the low 
energy limit, and up to first order in terms of the Gauss-Bonnet coupling a. In the limit a — > 
these expressions coincide with those extracted in [19]. However, in the presence of the Gauss- 
Bonnet term the implications are significant. Firstly, 4D holographic dark energy p\, apart from 
the usual squared holographic term, acquires a quartic correction. Secondly, matter density p is 
coupled with a holographic term oc R^ 2 , which is a result of p-p\ interaction of equation (3.14). 

Analytical solution of equations (3.19)-(3.21), namely finding H(a), then R^a), and finally 
PA( a ), is impossible. However, we are not interested in investigating the complete evolution but 
only in revealing the form of p\(a). Thus, we generalize Li's steps to construct a differential 
equation using Q\ as the unknown function [1]. 
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Firstly, we insert the usual variables: £Ia = jjfjrgi , &m = • Relation (3.19) then gives: 

ci 



Rh 



(3.22) 



up to C(a 2 ), with ci = and C2 = 12 ^ - Inserting this form in (3.21) and using the variable 
x = In a we obtain: 



00 dx _l 
Ha a 



+ ac 2 Jn A 'H] 



/Q^H 



(3.23) 



Similarly, using Qa, and Rh from (3.22), Friedmann equation (3.14) (with V±(a, PA5) given 
by (3.17)) up to 0{c?) writes: 



1 - n A = n M (1 + a 2c 3 n A H 2 ) 



(3.24) 



where C3 = 32-7r/3. In order to proceed forward we have to assume an explicit 0, m (o) dependence. 
In the interacting case at hand this should be different from the known ~ a -3 behavior of 
standard cosmology. However, in our model the p-pA interaction is downgraded by the extra- 
dimensional size as can be seen in (3.13) or equivalently in (3.17). Therefore, the deviation from 
conventional evolution will not be significant and we can use Qm = ^i^H^H^a^ 3 with 
and Hq the present values. Thus, we obtain: 

1 y/ay/l — TLa 



Ha 



tt M H 



i-ac 3 n A 

a 6 (I - Ua) 



(3.25) 



Finally, substituting this relation to (3.23) and taking derivative with respect to x, up to 0(a 2 ) 
we acquire the following differential equation: 



£l' A = Qi(Oa) + aQ 2 (flA,a) 



where 



Q 1 (£i A ) = nl(i-n A ) 



1 



+ 



and 



Q° H 2 \ 

Q 2 (n A ,a) = ^^{(c 2 -c 3 c 1 ] 



-5ni + Q 1 {n A ) (j^- 



2c 3 n 5 / 2 



(3.26) 



(3.27) 



(3.28) 



and the prime denotes the derivative with respect to x. Note that in the limit a — ► 0, differential 
equation (3.26) tends smoothly to that obtain by Li in [1], namely = Qi(0>a), an d can be 
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easily solved analytically. In the a ^ case of the present work such an exact solution is 
impossible. However, under the identification /3a (a) ~ a~ 3( - 1+WA \ we can extract the form of 
w\(z) at late times, i.e. at small z, with z = ^ — 1 and ao the value of a at present time (for 
simplicity we set ao = 1). We proceed as follows: 
Firstly, expanding In pa we obtain: 



In p A = lnp A |o + j,_ _ |o lna + - jn _ _", | (In a) 2 + ((In a) 3 ) 



(3.29) 



din a IU ' 2d(lna) 2 

where the derivatives are calculated at the present time ao = 1 [1]. Therefore, through p\(a) 
a — 3(i+«>a) we ma ke the identification: 



d In pa | . 1 d 2 lnp A | , n ,r 

a ma 2a(ma) z 



(3.30) 



Now, using Friedmann equation (3.24), and the expressions VL\ = 3 ^f^a and VLm 
n° M H^H- 2 a- 3 , we find: 



PA 



8vr a 3 (l - n A ) 



a 3 (l-fZ A ) 



(3.31) 



up to 0(a 2 ). Therefore, differentiating this relation with respect to lna = x, and using (3.26) 
for the calculation of the derivatives, we finally obtain the following u>a expression: 



wa(z) = w + w\z + a(w 2 + w 3 z), 



(3.32) 



where 



w 



3ci 



A- 



(3.33) 



Wl 



6ci 



1 + 



2,/0 A 



ci 



(3.34) 



i/' 2 



3ci 1 - ^9 



6ici + 26 2 6 3 ci 



i>i + b 2 b 3 ~ C1C362) 



(3.35) 



w 3 



O 
il A 



6c 2 1 



I -4(6i + 2& 2 6 3 )c 2 + ci(76i + 156 2 6 3 - 36 2 cic 3 )y / ^ A "+ 
+ (86 2 cic 3 - 661 - 8b 2 b 3 )n° A + ci[h - b 2 (3b 3 + c lC3 )}(n° A f/ 2 + 

+ 2[b 1 + 2b 2 (b 3 -c 1 c 3 )}(n° A ) 2 }. 



(3.36) 
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In the expressions above we have used the constants b\ = 2C3C4, 62 — c^/ci and 63 = C2 — C3C1, 
where C4 = ft^H^. Moreover, since ao = 1, we have replaced In a = — ln(l + z) ~ — z. Finally, 
0° is the present value of Oa- 

Relation (3.32) is the main result of this work and provides the Gauss-Bonnet correction to 
the corresponding result of [19]. Both investigations are formulated in the framework of bulk 
holographic dark energy. Therefore, although in the limit a — > 0, (3.32) coincides with Li's 
expression in [1], namely w\(z) = wo + w%z, the physical explanation in the present case comes 
through the 5D holographic consideration. This is the reason of the difference in constants 
between this work and [1]. 

From (3.32) it becomes obvious, that according to the value of c which is present in pas- 
relation (3.16), of C4 and of the Gauss-Bonnet coupling constant a, one can obtain a 4D holo- 
graphic dark energy behaving as phantom [50], quintessence or quintom [21, 22], i.e crossing the 
phantom divide j«a = — 1 [23, 43] during the evolution. Additionally, one can use observational 
results concerning dark energy evolution [24, 25] in order to estimate the bounds of the constant 
c of [1], i.e the bounds of c\ of the present work. In particular, observational data from type 
la supernovae give the best-fit value c\ = 0.21 within 1-a error range [51], while those from the 
X-ray gas mass fraction of galaxy clusters lead to c\ = 0.61 within 1-a [52]. Similarly, com- 
bining data from type la supernovae, cosmic microwave background radiation and large scale 
structure give the best-fit value c\ = 0.91 within 1-a [53], while combining data from type la 
supernovae, X-ray gas and baryon acoustic oscillation lead to c\ = 0.73 as a best-fit value within 
1-a [54]. Inserting this range of c\ values into our model one finds that wq < — 1 and w\ > 0, 
thus, within 1-a, he obtains a quintom-type holographic dark energy. Furthermore, W2 < 
while 1^3 > and therefore the Gauss-Bonnet contribution decreases the present value of w\, 
while it increases the growing rate of w/\_(z) with z, in comparison with the case where such a 
term is absent. However, the quantitative correction of the case will be very small, for 

reasonable C4 values. The reason is that, as we have mentioned, the p-p\ coupling, which arose 
naturally as a term V\{a, pas) P in (3.14), is downgraded by the extra-dimensional size as can 
be seen in (3.13) or equivalently in (3.17) (where we acquire a I? in the denominator). Thus, 
making the assumption that L5 is arbitrary large we downgrade the Gauss-Bonnet correction, 
too. It should be interesting to investigate the case where the bulk-size is smaller than the future 
event horizon, as in the two-brane model of [20, 55], but with the inclusion of a Gauss-Bonnet 
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term. The subject is under investigation. Finally, note that the role of the Gauss-Bonnet term 
on the w = —1 crossing has been investigated both in conventional 4D [32] and in braneworld 
frameworks [29, 31]. The novel feature of our work is the combined investigation of such a term 
with the bulk holographic dark energy. 

IV. DISCUSSION-CONCLUSIONS 

In this work we apply bulk holographic dark energy in a general braneworld model, with an 
induced gravity term and a perfect fluid on the brane, and a Gauss-Bonnet term in the bulk. 
Such a generalized bulk version of holographic dark energy is necessary if we desire to match 
the successes of brane cosmology in both theoretical and phenomenological-observational level, 
with the successful, simple, and inspired by first principles, notion of holographic dark energy 
in conventional 4D cosmology. In particular, as we showed in [19], the bulk space is the natural 
framework for the cosmological application, concerning dark energy, of holographic principle, 
since it is the maximally-dimensional subspace that determines the properties of quantum-field 
and gravitational theory, and the black hole formation. Subsequently, this bulk holographic 
dark energy will give rise to an effective 4D dark energy with "inherited" holographic nature, 
and this one will be present in the effective Friedmann equation. 

Taking the Gauss-Bonnet combination into account, a p-p\ coupling appears in the Fried- 
mann equation of the brane. We mention that this term arises naturally and is not a result 
of an inclusion by hand, which is the usual case of 4D interacting holographic dark energy in 
the literature [39, 47, 48]. This fact makes bulk holographic dark energy in the Gauss-Bonnet 
framework an interesting subject for further investigation. 

Examining the low-energy (late-time) evolution of the aforementioned model, we acquire 
the relation of w\(z) up to 0{a 2 ) and 0{z 2 ). In the limit a — > we re-obtain the results 
of [19] and those of conventional 4D calculations [1, 2, 3, 4, 5, 6], although in the 5D study 
the interpretation and explanation of these results is fundamentally different. In the presence 
of Gauss-Bonnet combination, and taking into account the constraints on the values of the 
constants by observational data, we find that the effective 4D holographic dark energy behaves 
as a quintom, i.e it crosses the phantom divide w\ = — 1 during the evolution. In particular, 
we observe that the presence of a non-zero a makes the current value of w\ smaller, while it 



12 



increases its growing rate with z, comparing to the a = case. However, the corresponding 
quantitative correction is very small due to the diminution of the p-p\ coupling by the arbitrary 
large extra-dimensional size. Yet, it should be interesting to investigate the case where the bulk 
size is smaller than the future event horizon. Then, the p-pA coupling would be significant and 
we would naturally acquire the advantages of interacting holographic dark energy, such as the 
coincidence problem solution, and the corresponding effects on wa(z). 
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